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Abstract 

A transport model based on the mean free path approach to describe pp colhsions is proposed. 
We assume that hadrons can be treated as bags of partons similarly to the MIT bag model. When 
the energy density in the collision is higher than a critical value, the bags break and partons 
are liberated. The partons expand and can make coalescence to form new hadrons. The results 
obtained compare very well with available data and some prediction for higher energies collisions 
are discussed. Based on the model we suggest that a QGP could already be formed in the pp 
collisions at high energies. 
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In recent works, we have proposed a kinetic approach to deal with a hot pion gas and a 
possible phase transition to the Quark Gluon Plasma (QGP)[l|. In our approach we treated 
the hadrons as bags of partons similar to the MIT bag model: if in a particle-particle collision 
the energy density is higher than a critical value obtained from the bag constant, new partons 
are created. Those partons can evolve in time and collide again with other particles. If the 
partons travel in a region whose density is below a critical value, they can coalesce and form 
new hadrons. The equation of state at zero barionic density was calculated and the effects 
of a mean field to reproduce Lattice QCD (LQCD) results were implemented as well[l|, Q]. 
Up to date LQCD results suggest that there is a cross over from a meson system to a QGP 
at a temperature of about 175 MeV Q]. 

It is the purpose of this work to extend the model to pp (and later to AA) collisions which 
are some of the most powerful tools to detect experimentally a transition to the QGP. Of 
course exact microscopic simulations for out of equilibrium-finite systems are out of reach 
at present. On the other hand transport approaches have been very useful in the past in 
describing many features of lower energies heavy ion collisions. Generalizations to relativistic 
energies of low energy heavy ion collisions ^, 0, (known as Boltzmann Uehling Uehlenbeck 
(BUU), Vlasov(VUU)/ Landau (LV)) have been proposed. The method we discuss in this 
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work is known as Boltzmann Nordheim Vlasov (BNV) approach at low energies 
is based on the concept of the mean free path approach 31. 

We can easily include the possibility of a QGP using the bag model [JJ, 19|, HOj. In fact, 
for each elementary hadron-hadron collision, we can calculate the local energy density and 
the pressure. If such quantities overcome the bag pressure and energy density, then Ug = rig 
quarks and antiquarks and Ug gluons are created. The number of quarks and gluons are 
calculated assuming local thermal equilibrium. In this way we can simulate a hadron gas 
and its transition to the QGP. In [1], we discussed the cases of Nj = 0,2,3, where Nf is 
the number of flavors plus a mean field. We also discussed the possibility that quarks could 
recombine and gluons decay into two quarks during the dynamical process. We saw that 
in order to have a reasonable description of hadronization the local quark (gluon) density 
must not exceed a critical density calculated from the MIT bag model. If the local density 
is larger than such a value (which is equivalent to having a high temperature or energy 
density) the quarks cannot recombine to form hadrons (or the gluons decay into qq pair). 
From a comparison of our results to LQCD we realized that in the QGP phase our e/T^ is 
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higher than the value suggested by lattice calculation and indeed it approaches the Stefan- 
Boltzmann limit as it should be. Knowing that there is Debye color screening in the QGP 
phase we calculated the corresponding mean field and adjusted the parameters to reproduce 
lattice data l| . 

In order to use the same approach for out of equilibrium situations we have to deal first 
with energetic pp (or ee) collisions. This is the purpose of our paper and we will show 
in the following that keeping the main ingredients of the model but assuming a suitable 
momentum distribution of the partons after they are created we can reproduce reasonably 
well some data on pp collisions. Of course, the EOS which can be obtained from the model 
remains unchanged by these modifications. Thus in our approach we can describe with some 
accuracy, both elementary collisions as well as collective properties of the systems. Quantum 
statistics (i.e. Pauh and Bose statistics) are included similarly to for Bose and [121] for 
Fermi statistics l|]. 

The number of particles that are produced in each elementary collision increases with 
increasing beam energy in agreement to data. We make some prediction for pp collisions 
at 5 TeV which will be soon studied experimentally at CERN. The model assumes that 
hadrons can be produced from partons coalescence similar to the results of low energy 
heavy ion collisions to study a liquid-gas phase transition We analyze our results in 
terms of the Fisher model for a liquid-gas phase transition and surprisingly we find a power 
law mass distribution with an exponent r = 2.3. This is an expected value from Fisher's 
approach[lJ]. For comparison, we have collected some data obtained in Au+Au collisions 
at RHIC and we found a similar mass distribution. This calls for a critical behavior of the 
system. We suggest that this could be due to the rapid expansion of the system after the 
partons have been created. If we push the analogy to the liquid-gas phase transition, we 
can assume that the vapor is made of partons and maybe pions, while the other hadrons 
make the liquid part. After the QGP is created, the system expands rapidly and since the 
partons cannot remain isolated they coalesce. This results in a rapidly falling (power law) 
mass distribution. 



I. FORMALISM 



The mean free path method discussed above has been studied in detail at low energies 
and it has been shown to solve the Boltzmann equation in the cases where an analytical 
solution is known [3]. We have generalized the approach to keep into account relativistic 
effects. The particles move on straight lines during collisions since we have not implemented 
any field in this exploratory study. For short we name the method proposed as Relativis- 
tic Boltzman equation (ReB). In order to test our approach, we have discussed also some 
simple cases where analytical solutions are known to verify the sensitivity of our numerical 
approximations [l, 7, 8|. 

In the present calculations we simulate pp collisions at a given energy and impact pa- 
rameter. For each energy, events are obtained changing the initial impact parameter. 




FIG. 1: Schematic view of the pp collision process, b is the impact parameter and r the relative 
distance between the two protons. 



We recall how we include partons in our approach. First, for a massless quark and gluon 
plasma in equilibrium the following relations hold for the pressure P, quark (antiquark and 
gluon) density and energy density e versus temperature T 
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where e = 3P, gtot = 16 + y^/' number of flavors. In the MIT bag model 
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quarks and gluons are confined in the bag if tlie pressure is less tlian tlie critical pressure B 
tliat tlie bag can sustain. Thus from the previous equation we can assume that the quarks 
and gluons are liberated in a collision if the energy density is larger then 3B. This gives 
a critical energy density ec = 3B = 1.69GeV fm using 

^1/4 ^ o.242Gel^|9|, [lo| and for 
three flavors. For each h-h elementary collision we know the energy of the collision and 
the interaction volume from the distance between the colliding particles in their center of 
mass system. This distance is taken as the radius of a sphere enclosing the two colliding 
particles and subsequently the newly formed partons. Thus we can calculate the number of 
quarks, antiquarks and gluons as a function of the energy density liberated in the collision 
inverting equation (1). These relations can be easily generalized to the cases of finite quark 
masses. In this work we will use the values of 5,10 and 160 MeV for u,d and s-quark masses 
respectively. We stress that these relations are strictly valid in thermal equilibrium but we 
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FIG. 2: Analytical estimate of the number of charged particles produced in pp collisions as function 
of energy. Data are given by the full squares, calculations are given by the open circles (without 
probability corrections) and open squares. 

are extending them here for non equilibrium cases, i.e. the relevant control parameter is the 
energy density and not the temperature. Of course, this is a strong ansatz and it could be 
justified if the system is very chaotic, thus many particles should be created for each collision, 
and also after averaging over many ensembles to have some statistical meaning. A crucial 



point of this approach is the momentum distribution just after the partons are created. For 
particles in a box the distribution did not really matter much, since, because of elastic and 
inelastic collisions equilibrium is finally reached. This is not true in pp collisions, in fact we 
know experimentally that in such cases equilibrium is not reached. Thus we assume that 
the partons created have originally a momentum given by a fraction of the initial z-beam 
momentum. We assume (as we did in that the partons can collide successively elastically 
and inelastically and we fit the corresponding cross sections to reproduce the pseudorapidity 
experimental data. 

If in a collision between a parton and a hadron the local energy density is larger than 
the critical value, new partons are liberated from the hadron similarly to the mechanism 
discussed above. After being created, the partons expand and coalesce into hadrons if the 
local density become smaller than a critical value obtained from the MIT bag model l| . 
In the model, coalescence occurs through resonances which decay later on[15]. We have 
included known resonances and their decay using parametrizations of the data up to 2 GeV. 

Given the essential ingredients of the model, the crucial point now is to reproduce for 
instance the number of produced particles as a function of energy. We can estimate analyt- 
ically the number of produced particles making some simplifications. Let us consider two 
colliding protons along the z-axis and at impact parameter b, assume the particles move 
on a straight line until they collide. We can estimate a critical volume Vc when the energy 
density is equal to the critical value for breaking the bags. The situation is depicted in figure 
(1), and we easily get: 



zq is the coordinate where, for a given b, the volume enclosing the two colliding particles 
is equal to the critical one, see fig.l. The average (over impact parameters) number of 
produced quarks is given by: 



and the maximum impact parameter bmax can be obtained from equation (2), zq = 0. 
The last term in the integration gives the probability that particles collide at a given point z: 
n(6) = 1 — e^o^^/-^ , and the mean free path A = We can further simplify the integration 




(2) 





by assuming that the particles colhde after travehng a distance zq, i.e. z=0. The estimate 
is compared to experimental data from [l^ in figure 2 and we can see that we obtain rather 
surprisingly a reasonable agreement. In particular, the agreement is improved if one takes 
into account the collision probability in eq.(2) (open squares in fig. 2). 
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FIG. 3: Cross section versus energy from the analytical estimate, eq(5). The data are given 
by the full squares, open squares are our results while the crosses are the input values in the 
probability,eq.(4) . 

In our simulation a collision between particles is probabilistic, different from other 
approaches 

mm 

where a geometrical assumption is made for collisions (black sphere). 
This implies that collisions (both elastically and inelastically), can occur as well for impact 
parameters larger than ^cr/n). Thus in order to know what is the effective cross section im- 
plemented we need to integrate over impact parameters. For instance, for inelastic collisions 
we get: 

f^max 

ai = 27i bU{b)db (4) 
Jo 

Notice that this equation for a contains a cross section also in the probability, in some 



7 



sense we are transforming a sharp cutoff probability into a diffuse one. Recall that the 
inelastic process depends from the reaching of a critical energy density in our model, thus 
from bmax, and this process will result in a certain dependence on energy of the cross section. 
Our strategy is to fit the input cross section in the probability in order to reproduce the 
experimental values. This is shown in figure (3) where we plot the experimental cross sections 
(full squares) versus energy 10(] . We reproduce rather well the data (open squares) when we 
choose an input cross section given by the crosses. We see that we are able to reproduce 
the results especially at high energies. The input cross section is one of the free parameters 
of our model together with the critical energy density and the quarks masses. Once we 
have fixed such values we obtain the number of charged particles shown in figure(2). In 
the numerical simulation that we discuss below, we do not make any assumption and the 
trajectory is followed in time without simplifications. It is also included the possibility that 
partons are created from secondary collisions during the expansion phase. For our purposes, 
the result found here gives a strong hint that the assumption we make are rather reasonable. 
However, we expected this result since we know that statistical models work well already in 
pp and ee collisions as claimed for instance in jl6| . 

II. NUMERICAL SIMULATIONS 

Numerically the collisions are followed in time for given initial conditions and beam 
energiesjl, 7|. We have performed simulations from tens of GeV initial energy to several TeV 
which should be soon experimentally available at CERN. The average number of produced 
charged particles and the corresponding cross sections are calculated by generating many 
events for each energy and for many impact parameters b. The number of events is generated 
proportionally to the impact parameter b as usual. eq.(3). In fi. (4) we plot the obtained 



cross section as function of y (s). A parametrization of the data [lO| is given by the full dots 
and our numerical results by the open dots. The experimental points at the highest energy 



were obtained in air shower experiments 17 1. 



The average number of produced charged particles can be calculated similarly to the 
previous section and the results are given in fig. 5. 

Because of the probabilistic nature of our approach the number of particles produced 
depend strongly on the impact parameter b. In figure (6) we plot the number of produced 
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FIG. 4: Cross section versus energy. The data is given by the fuh dots, open dots are our numerical 
results. 

charged particles versus b/bmax at different beam energies. Since the maximum impact 
parameter increases with energy, see eq.(2), we reproduce the well known Blacker, Edger 
and Larger (BEL) effect 3 : for higher energies the elementary cross section increases thus 
the number of involved impact parameters increases as well. 

A more sensitive quantity to obtain is for instance the pseudorapidity distribution of 
charged particles. In our approach after the partons are created they are given an initial 
random momentum along the z-axis (the beam axis) and proportional to the initial momenta 
of the colliding hadrons. Total energy and momenta conservation are enforced. In the 
following time steps, the partons expand and they could collide with other particles getting 
in this way some momentum in the transverse direction to the beam axis. Also diquarks can 
be formed with a given probability. The diquarks can subsequently collide again and form a 
hadron. We fit the parton-particle and diquark-particle cross sections in order to reproduce 
the experimental data on pseudorapidity distributions [lo|, [l^. In fig. 7 the distribution 
is plotted at -y/(s) = 2006*6^^ for different parton-particle cross sections. Notice that the 
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FIG. 5: Number of produced charged particles versus energy. The data is given by the fuh symbols. 

transverse momentum is obtained not only from these collisions but also from the coalescence 
and further decay of resonances later on when the density is smaller than a critical value as 
obtained from the MIT bag model. However, we notice that a cross section of few mb works 
rather well. In figure 8 we plot the distributions for different initial energies. The agreement 
to data is fair. 

In fig. (9) we have performed calculations at energies larger than 1 TeV which is our 
prediction for some future CERN experiments. 

III. MASS DISTRIBUTIONS 

The model we have proposed makes use of a coalescence mechanism to produce the 
new hadrons. The physics we have in mind is that after the pp collision, partons are 
created in a very dense region. The partons expand and eventually collide thus getting some 
transverse momentum. When the local density decreases below a certain value, partons 
coalesce into resonances. Those hadrons expand further and eventually decay or collide 
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FIG. 6: Number of produced charged particles versus b at different energies. The dashed hues give 
the average number (over impact parameters) of produced charged particles. 

with other particles. This mechanism is reminiscent of droplets formation at lower energies 
for a liquid to gas phase transition 13]. In the same sense we could compare the partons 
to the vapor and the hadrons to the liquid. In the expansion the system crosses from one 
phase to the other. In the crossing we expect fluctuations to be large which could be reveled 
from the fluctuations of momenta distributions of hadrons 20|| . A strong signature which 
appears also in non equilibrium cases and for a reduced number of particles is a power 



law in the mass distributions 



19|]. In a pp collisions we have seen that the multiplicity 



of produced particles could be large, thus we could expect in the model that a self-similar 
mass distribution appears. In fig. (10) we plot the numerical mass distributions Y{> m) 
which is the number of fragments obtained with mass greater than m[2]|. We have adopted 
this definition rather than the commonly used yield to avoid numerical fluctuations in the 
simulations. We can fit such a distribution with the expression: Y{> m) = cost x m^'^^^. 
From the fit we obtain a value of r = 2.3 which is very reminiscent of the Fisher law for 
a gas to liquid phase transition. Thus this result hints for a second order phase transition. 
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FIG. 7: Pseudorapidity distribution at ^/s = 200Gey, the data points are given by the full dots 18l |. 



Theoretical calculations are given by the open symbols for different parton-particle cross sections. 

As we see from the figure, the power law is rather strong and independent on the beam 
energy, the rapid falloff at around 2 GeV is simply due to the inclusion in the calculations 
of resonances up to 2 GeV and to the low statistics. 

We expect that this result is independent on the way the plasma was formed, thus should 
show even better in nucleus- nucleus collisions. We have collected some experimental data 
for Au+Au collisions at RHIC 2J] which are displayed in fig. (11). In this case since the data 
has enough statistics we have plotted directly the measured yield. Also notice that now the 
power law is extended to masses larger than 2 GeV. The bumps and dips in the figure are 
due to detected resonances of different widths. This is somehow similar to the corresponding 
oscillations in the mass distributions at lower energy for liquid-gas phase transition due to 
shell effectsQ 
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FIG. 8: Pseudorapidity distributions at different beam energies. Data are given by the full 
symbols 

IV. CONCLUSIONS 

In this work we have applied a recently introduced transport approach to study proton- 
proton collisions. This is a necessary step if we want to simulate nucleus-nucleus collisions. 
The model includes resonance formation and their decays. The possibility of a QGP is in- 
cluded based on the MIT bag model as well as quantum statistics. We have seen that we 
can reasonably reproduce available data on pp collisions and also be able to make prediction 
at LHC energies. The main assumption, we believe, is inspired by thermal models which 
we have generalized to non-equilibrium situations. The success in the data reproduction 
suggests that the system is chaotic enough to justify our assumption. Pseudorapidity dis- 
tributions can be reproduced assuming some suitable parton-parton collision cross section. 
Hadronization occurs via partons coalescence which is similar to coalescence of vapor into 
drops. We have explored this analogy by analyzing physical quantities which are relevant 
in phase transitions such as mass distributions and found a power law with an exponent 
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FIG. 9: Pseudorapidity distributions at 1.8,5 and 10 TeV. 

compatible with a Fisher's law analysis. This is seen not only in our calculations but also 
in Au+Au experimental data. This signature is at the moment not sufficient to pin down 
the phase transition and should be considered as a curious result. In following works we 
will present a detailed analysis for a critical behavior and we hope to see soon corresponding 
data to con&m or reject our scenario. 
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formed). We also thank V.Greco for stimulating discussions and comments. 
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